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Abstract 

The precise determination of the Z boson parameters from the measure- 
ments at the Z resonance by the four collaborations Aleph, Delphi, L3 and 
Opal in e + e~ collisions at the large electron positron collider Lep at CERN is 
a landmark for precision tests of the electroweak theory. The four experiments 
measured quantities which were used to extract the mass and width of the 
Z boson, the hadronic cross-section at the pole of the resonance, the ratio of 
hadronic and leptonic decay widths, and the leptonic forward-backward asym- 
metries at the pole. The combination procedure based on these parameters is 
presented in this paper. 

(to be published as part of a review in Physics Reports) 
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1 Introduction 



Between the years 1989 and 1995 the e + e~ collider Lep at CERN provided interactions 
at centre-of-mass energies, yfs, ranging from 88 to 95 GeV, i. e. around the mass of the 
Z boson (Lep I phase). An important aspect of physics at Lep concerns the analysis 
of fermion-pair production in e + e~ collisions The four Lep experiments Aleph ||, 
Delphi Jt|, L3 || and Opal || analysed, in particular, hadron (quark-pair) production 
and the pair production of charged leptons, £ — e, //, r. 

At various centre-of-mass energies, total cross-sections are measured for all processes, 
while forward-backward asymmetries are measured in lepton-pair production. These mea- 
surements ("realistic observables") allow the determination of various properties of the Z 
boson such as its mass, total and partial decay widths, and coupling constants to fermions 
("pseudo-observables" ) . For the extraction of the pseudo-observables, the experiments per- 
form model- independent fits to their measured realistic observables 

To obtain the best possible precision the results of the four Lep experiments have 
to be averaged. This paper describes the combination procedure adopted by the Lep 
electroweak working group. Performing an average over the realistic observables constitutes 
an extremely complicated task, as it involves hundreds of measurements, each with specific 
phase space definitions and experimental errors which are correlated among different centre- 
of-mass energies and data taking periods and also among the Lep experiments. Therefore, 
the combination of the experimental results is performed on the basis of the four sets 
of pseudo-observables. As will be shown here, this is possible without significant loss of 
precision. 

Additional fits to the experimental data, usually not contained in the individual publi- 
cations quoted above, were provided by the experiments and are documented in this paper. 
A large effort is devoted to the treatment of systematic errors and their correlation among 
the experiments, such that the combination procedure yields an optimal estimator for the 
averages. 

This paper is organised as follows: Section^ summarises the information about the 
Lep I running relevant for this paper. Section^ presents a brief introduction to the pseu- 
do-observables used to parametrise the realistic observables around the Z resonance. The 
individual experimental results are presented in Section |. Sources of correlated systematic 
errors between experiments and their effects on the pseudo-observables are discussed in 
Section^. In Section^, studies of various methods for combining the results are presented. 
The resulting pseudo-observables are then considered in the framework of a specific model, 
the minimal Standard Model (SM), and are compared with direct SM fits to the realistic 
observables. The conclusions are summarised in Section |^. 

2 Z resonance scans at LEP I 

Running of Lep in the years from 1989 to 1995 was dedicated to precision studies of the 
Z boson parameters. Electron-positron collisions were provided at several well-determined 
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centre-of-mass energies around the Z resonance, with steadily improving performance. 
The set of measurements collected by the experiments consists of the hadronic and lep- 
tonic cross-sections and the leptonic forward-backward asymmetries around seven points 
in centre-of-mass energy, over six years of running at Lep I. In addition, changes of experi- 
mental conditions, such as the inclusion of new detector components, made it necessary to 
subdivide the data samples even further. The full Lep I data set consists of about 4 x 200 
individual cross-section and asymmetry measurements. From these each experiment has 
extracted a set of parameters describing the cross-section around the Z resonance, which 
include the mass, mz, and width, Tz, of the Z and the total pole cross-section for qq pro- 
duction, These parameters are discussed in detail in Section QED corrections from 
initial-state photon radiation are large around the Z resonance due to the rapid variation of 
the cross-sections with centre-of-mass energy. For illustration, Figure [I] shows the average 
over the hadronic cross-section measurements by the four experiments, together with the 
fitted line shape curve before and after unfolding photon radiation. The cross-section is 
dominated by on-shell Z production, although photon exchange and 7-Z interference con- 
tributions are not negligible. The measurements are sensitive to higher-order electroweak 
corrections. These modify the tree-level couplings of the Z to fermions, and are quantified 
in terms of electroweak form factors. 

Much effort was dedicated to the determination of the energy of the colliding beams, 
which reached a precision of about 20 x 10~ 6 on the absolute energy scale. This level of ac- 
curacy is vital for the measurements of the mass and width of the Z. All the experiments 
replaced their first-generation luminosity detectors, which had systematic uncertainties 
around the percent level, by high-precision devices capable of pushing systematic errors 
on the acceptance of small-angle Bhabha scattering events below one per-mil. As a conse- 
quence of improvements of the accelerator and of the experiments during Lep I running, 
statistical and systematic errors are much smaller for the last three years of data taking, 
which hence dominate the precision achieved on the Z parameters. 

2.1 Event selection and statistics 

During the summer of 1989 the first Z bosons were produced at Lep and observed by 
the four experiments. Since then the operation of the machine and its performance were 
steadily improved. At the end of data taking around the Z resonance in autumn 1995 the 
peak luminosity had reached nearly twice its design value. Table|l| summarises the data 
taking periods, the approximate centre-of-mass energies and the integrated luminosities 
delivered. 

The data collected in 1989 constitute only a very small subset of the total statistics and 
are of lower quality, and therefore are not used here. In the years 1990 and 1991 "energy 
scans" were performed at seven different centre-of-mass energies around the peak of the 
Z resonance, placed about one GeV apart. In 1992 and 1994 there were high-statistics 
runs at the peak energy only. In 1993 and 1995 data taking took place at three energy 
points, about 1.8 GeV below and above the peak and at the peak. In particular the off-peak 
energies were carefully calibrated employing the technique of resonant depolarisation of the 
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Figure 1: Average over measurements of the hadronic cross-sections by the four exper- 
iments, as a function of centre- of -mass energy. The dashed curve shows the QED de- 
convolved cross-section, which defines the Z parameters described in the text. 



transversely polarised beams ||10Hl3|]. 



The accumulated event statistics amount to about 17 million Z decays recorded by the 
four experiments. A detailed breakdown is given in Table 0. 

As an example, the measurements of the hadronic cross-section at the three principal 
energy points are shown in Figure 0. Because the hadron statistics are almost ten times 
larger than the lepton statistics, these measurements dominate the determination of the 
mass and the width of the Z. Detailed descriptions of the individual experimental anal- 
yses can be found in |[]-§J. They all rely on excellent separation of the final states, qq, 
e + e~, and r + r~, accompanied by high selection efficiencies. The total cross-section, 

o~ tot , is determined from the number of selected events in a final state, N, the number of 
background events, N hg , the selection efficiency, e, and the integrated luminosity, C, as 
a tot = (N-N hg )/(eC) . 
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Figure 2: Measurements by the four experiments of the hadronic cross-sections around the three principal energies. The 
vertical error bars show the statistical errors only. The open symbols represent the early measurements with typically 
much larger systematic errors than the later ones, shown as full symbols. Typical experimental systematic errors on the 
determination of the luminosity are indicated in the legend; these are almost fully correlated within each experiment, but 
uncorrelated among the experiments. The horizontal error bars show the uncertainties in Lep centre- of-mass energy, 
where the errors for the period 1993-1995 are smaller than the symbol size in some cases. The bands represent the result 
of the model-independent fit to all data, including the two most important common theoretical errors from initial-state 
photon radiation and from the calculations of the small-angle Bhabha cross-section. 
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[pb- 1 ] 

Ll J 


1989 


[88.2,94.2] 


1.7 


1990 


[88.2,94.2] 


8.6 


1991 


[88.5,93.7] 


18.9 


1992 


91.3 


28.6 


1993 


89.4,91.2,93.0 


40.0 


1994 


91.2 


64.5 


1995 


89.4,91.3,93.0 


39.8 



Table 1: Approximate centre- of-mass energies and integrated luminosities delivered by 
Lep, per experiment. In 1990 and 1991, a total of about 7 pb^ 1 was taken at off-peak 
energies, and 20 pb^ 1 per year in 1993 and in 1995. The total luminosity used by the 
experiments in the analyses was smaller by 10-15% due to data taking inefficiencies. 
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year 


A 


D 


L 





all 


year 


A 


D 


L 





all 


'90/91 


433 


357 


416 


454 


1660 


'90/91 


53 


36 


39 


58 


186 


'92 


633 


697 


678 


733 


2741 


'92 


77 


70 


59 


88 


294 


'93 


630 


682 


646 


649 


2607 


'93 


78 


75 


64 


79 


296 


'94 


1640 


1310 


1359 


1601 


5910 


'94 


202 


137 


127 


191 


657 


'95 


735 


659 


526 


659 


2579 


'95 


90 


66 


54 


81 


291 


total 


4071 


3705 


3625 


4096 


15497 


total 


500 


384 


343 


497 


1724 



Table 2: The qq and £ + £ event statistics, in units of 10 3 7 used for the analysis of the 
Z line shape and lepton forward-backward asymmetries by the experiments Aleph (A), 
Delphi (D), L3 (L) and Opal (0). 

The luminosity of the beams is determined by normalisation to the theoretical cross- 
section for the process of small-angle Bhabha scattering, which is dominated by photon 
exchange in the t channel. Thus, the integrated luminosity is given by the number of 
observed small-angle Bhabha events and the calculated cross-section for this process within 
a given experimental acceptance. This requires the detection of back-to-back electrons 
and positrons close to the beam direction. Their positions and energies are precisely 
measured by forward calorimeters placed at small angles with respect to the beam line 
(typically 30 mrad< 6 < 50 mrad). The main experimental systematic error arises from 
the definition of the geometrical acceptance for this process. Since the angular distribution 
is steeply falling with increasing scattering angle (oc 6~ 3 ), the precise definition of the inner 
radius of the acceptance region is most critical. 

The forward-backward asymmetry, Ar>, is defined by the numbers of events in which the 
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final state lepton goes forward or backward with respect to the direction of the incoming 
electron, Nf and N h , respectively: Ar> = (Nf— N h )/ (Nf+N h ) . In practice, Ar, is determined 
from a fit to the differential cross-section of the form da jd cos 9 oc 1 + cos 2 9 + 8/3 Ar, cos 9, 
where 9 is the angle between the direction of the final state fermion and that of the incoming 
electron. This procedure makes better use of the available information and hence leads 
to slightly smaller statistical errors. The electron final state is special due to t-channel 
diagrams, as is discussed in more detail in Section |572| . The forward-backward asymmetries 
do not require any normalisation, but rely on precise measurements of the production 
angles of the final state fermions. Forward-backward asymmetries in qq final states are not 
considered here, because these require either dedicated techniques for the tagging of quark 
flavours or a special method to extract the inclusive quark forward-backward asymmetry 
from the natural mixture of quark flavours in hadronic events. 

In general, the systematic errors arising from the selection procedures are small and so 
the accumulated statistics can be fully exploited. Furthermore, the purely experimental 
errors arising from the limited understanding of detector acceptances are uncorrelated 
among the experiments. Errors arising from limitations in theoretical precision, such as the 
calculation of the small-angle Bhabha cross-section, the t-channel contribution in the e + e~ 
final state or pure QED corrections to the cross-section, are common to all experiments. 
They are discussed in detail in Section[|. 

An overview of the experimental systematic errors is given in Table|3|. The systematic 
error on the luminosity is common to cross-section measurements of all final states, but 
does not affect the measurements of Ar,. 



2.2 Energy calibration 

Precise knowledge of the centre-of-mass energy is essential for the determination of the 
mass and width of the Z resonance. The key features of the energy calibration procedure 
were the technique of resonant depolarisation and the careful monitoring of all relevant 
machine parameters ||13|| . The latter is important because beam energy calibrations with 
resonant depolarisation were possible only outside normal data taking, usually at the end 
of data taking in a particular fill of the accelerator, with fills typically lasting approxi- 
mately 10 h. About 40% of the recorded off-peak luminosity was calibrated in this way 
in the 1993 scan and about 70% in the 1995 scan. For each experiment a value of the 
beam energy was provided every 15 minutes. These values were evaluated from the time 
evolutions of the relevant machine parameters. This required a model which took into 
account the fields in the Lep dipoles and in the corrector magnets, beam orbit positions, 
collision offsets at the interaction points and parameters of the radio-frequency system. In 
addition, environmental effects from leakage currents produced by trains in the Geneva 
area and the gravitational forces of the Moon and the Sun leading to small deformations 
of the accelerator geometry had to be considered. Errors on the centre-of-mass energy 
are largely dominated by the uncertainties in this model. The energy errors vary slightly 
among the interaction points, largely due to different configurations of the radio frequency 
cavities. The energy errors for different experiments and periods of data taking have large 
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Aleph 






Delphi 






'93 


'94 


'95 


'93 


'94 


'95 




0.067% 


0.073% 


0.080% 


0.24% 


0.09% 


0.09% 


0"had 


0.069% 


0.072% 


0.073% 


0.10% 


0.11% 


0.10% 


0~ e 


0.15% 


0.13% 


0.15% 


0.46% 


0.52% 


0.52% 


°> 


0.11% 


0.09% 


0.11% 


0.28% 


0.26% 


0.28% 


a T 


0.26% 


0.18% 


0.25% 


0.60% 


0.60% 


0.60% 


ID 


0.0006 


0.0006 


0.0006 


0.0026 


0.0021 


0.0020 


/In, 
ID 


0.0005 


0.0005 


0.0005 


0.0009 


0.0005 


0.0010 


ID 


0.0009 


0.0007 


0.0009 


0.0020 


0.0020 


0.0020 






L3 






Opal 






'93 


'94 


'95 


'93 


'94 


'95 


£exp. 


0.086% 


0.064% 


0.068% 


0.033% 


0.033% 


0.034% 


0"had 


0.042% 


0.041% 


0.042% 


0.073% 


0.073% 


0.085% 




0.24% 


0.17% 


0.28% 


0.17% 


0.14% 


0.16% 


°» 


0.32% 


0.31% 


0.40% 


0.16% 


0.10% 


0.12% 


Or 


0.68% 


0.65% 


0.76% 


0.49% 


0.42% 


0.48% 




0.0025 


0.0025 


0.0025 


0.001 


0.001 


0.001 




0.0008 


0.0008 


0.0015 


0.0007 


0.0004 


0.0009 




0.0032 


0.0032 


0.0032 


0.0012 


0.0012 


0.0012 



Table 3: Experimental systematic errors for the analysis of the Z line shape and lepton 
forward-backward asymmetries at the Z peak. None of the common errors discussed in 
Section^ are included here. 

common parts, and therefore the use of a full correlation matrix is necessary. Assuming 
that all experiments contribute with the same weight allows all the Lep energy errors to 
be conveniently summarised in a single error matrix, common to all interaction points, as 
given in Table 

The energy of individual beam particles is usually not at the mean value considered 
above, but oscillates around the mean beam energy. Therefore observables are not mea- 
sured at a sharp energy, E® m , but instead their values are averaged over a range in energies 
E® m ± 5E cm . With the assumption of a Gaussian shape of the energy distribution, the 
total cross-sections receive a correction proportional to SE^ m and the second derivative of 
a(E cm ) w.r.t. E cm . At Lep I, typical values of the centre-of-mass energy spread are around 
50 MeV. The effects of the correction lead to an increase of the cross-section at the peak of 
the Z resonance by 0.16% and a decrease of the width by about 5 MeV. The uncertainties 
on the energy spread, around ±1 MeV in 1993-1995, constitute a negligible source of error 
common to all experiments. 

In addition to the natural energy spread, changes in the mean beam energy due to 
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'93 


'93 


'94 


'95 

■ / 'J 


'95 

i J 'J 


'95 

i J 'J 




p-2 


P 


p+2 


P 


p-2 


P 


p+2 


'93 p-2 


3.42 














'93 p 


2.76 


6.69 












'93 p+2 


2.59 


2.64 


2.95 










'94 p 


2.25 


2.38 


2.16 


3.62 








'95 p-2 


1.29 


1.14 


1.23 


1.23 


1.78 






'95 p 


1.19 


1.20 


1.25 


1.30 


1.24 


5.39 




'95 p+2 


1.20 


1.15 


1.33 


1.24 


1.22 


1.34 


1.68 



Table 4: Signed square root of covariance matrix elements, (V E ), in MeV ; from the deter- 
mination of the centre- of-mass energies for the scan points in 1993-1995 /|73|/. Elements 
above the diagonal are left blank for simplicity. The errors for the earlier years may be 
found in References\T(\, [ill. 



changes of machine parameters have a similar effect. Periods of data taking with a very 
similar centre-of-mass energy were combined into a single energy point in the experimental 
analyses by performing a luminosity-weighted average. The additional energy spread re- 
sulting from this was only around 10 MeV, which adds in quadrature to the natural energy 
spread of the accelerator. 

3 Parametrisation of the differential cross-section 

The differential cross-section for fermion pair production around the Z resonance consists of 
three s-channel contributions: from Z exchange, photon exchange and from the interference 
between photon and Z diagrams, a(s) = a z + a 1 + a mt . This can be cast into a Born-type 
structure with complex-valued, s-dependent form factors describing the couplings of the 
Z and the photon to fermions. In the Z pole approximation, valid for s ~ m|, these are 
taken to be constants. Neglecting initial and final state photon radiation, final state gluon 
radiation and fermion masses, the electroweak kernel cross-section can thus be written as 



2 s 1 da cw 
7r N { c d cos 6 



e e 



ff) 



a(mz)Q { (1 + cos 2 



7 



8Re{a*(m z )Q f X (s) [G e v G l v {l + cos 2 6) + 1Q\Q\ cosfl]} 
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Z interference 



+16 |x(*)| 2 [ (|^| 2 + m 2 ){\Q'v\ 2 + l^| 2 )(l + cos 2 ■ 
~ +8Re{^£nRe{6v£i*}cos0 
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with 




X(s) 



8irV2 s-ml + isr z /. 



Here a(mz) is the electromagnetic coupling constant at the scale of the Z mass, Gf is the 
Fermi constant, Q f is the charge of the final state fermion, and the colour factor N* is one 
for leptons (f=e, /i,r) and three for quarks (f=u, d, s, c, b). The effective vector and axial 
vector couplings of fermions to the Z are denoted by Q v and Q\. x( s ) is the propagator 
term characterized by a Breit-Wigner denominator with an s-dependent width. 

In Bhabha final states, e + e~ — > e + e~, the t-channel diagrams also contribute to the 
cross-sections, with a dominant contribution at large cos#, i.e. in the forward direction. 
This contribution, as well as its interference with the s-channel, add to the pure s-channel 
cross-section for e + e~ — > e + e~. 

The 1 + cos 2 9 terms in the above formula contribute to the total cross-section, whereas 
the terms multiplying cos 9 contribute only to the forward-backward asymmetries for an ex- 
perimental acceptance symmetric in cos 9. The total cross-section is completely dominated 
by Z exchange. The 7 — Z interference determines the energy dependence of the forward- 
backward asymmetries and dominates at off-peak energies, but the leading contribution 
from the real parts of the couplings vanishes at y/s = m Z - 

The inclusion of higher-order electroweak corrections is absorbed in Qa and Qy with 
small imaginary parts arising from electroweak form factors. The experimental measure- 
ments do not allow a simultaneous extraction of the real and imaginary parts, and therefore 
the effective couplings to be determined are defined as qa = Re(^) and gy = Re(Qy). 
The imaginary parts of Qa and Qy are explicitly accounted for in the fitting codes by 
setting them to their SM expectations. The effects of box diagrams are also taken into 
account at this level. 

It is worth noting that the definitions of the mass and width with an s-dependent width 
term in the Breit-Wigner denominator are suggested by phase-space and the structure of 
the electroweak radiative corrections within the SM. They differ from other commonly used 
definitions, e. g. the real part of the pole position in the energy-squared plane, where the 
propagator term takes the functional form x( s ) ^ s /( s ~ ™z 2 + i^z^z) ■ This gives an 
identical line shape if mz and Tz are related to mz and Tz by the multiplicative factor 



Photon radiation from the initial and final states, and their interference, is conveniently 
treated by convoluting the electroweak kernel cross-section, cr ew (s), with a QED radiator, 



The difference between the forward and backward cross-sections entering into the determi- 
nation of the forward-backward asymmetries, op — &b, is treated in the same way using a 
radiator function #qed- 

At the peak the QED de-convoluted total cross-section is 36 % larger than the measured 
one, and the peak position is shifted downwards by about 100 MeV. The estimated precision 



Vv/i + r 2 z /mi ■ 
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of this important correction is discussed in Section p.4.1| . 

The partial Z decay widths are defined inclusively, i. e. they contain QED and QCD final 
state corrections and the contribution from the imaginary parts of the effective couplings, 

Ffi = N* ^ F / _ Z (\Ga\ 2 Ra + \Gv\ 2 Rv) + ^ew/QCD- 



6^7T 

The radiator factors R v and R A take into account final state QED and QCD corrections as 
well as non-zero fermion masses; A cw /qcd accounts for non-factorizable electroweak/QCD 
corrections. The inclusion of the imaginary parts of the couplings in the definition of the 
leptonic width, leads to changes of 0.15 per-mil corresponding to only 15% of the 
LEP-combined experimental error on T^. 

The total cross-section arising from the cos ^-symmetric Z production term can also be 
written in terms of the partial decay widths into the initial and final states, r ee and Tg, 



a s = a. 



„r2 

peak a 1 Z 



2r2 , „ ' 



• ,, peak 1 j 12 7T T^Fff 

with < = — — a s and a s = — T — . 

1 + Cqed ra z 1 z 

The term 1/(1 + £qed) removes the final state QED correction included in the definition 

of r ee . 

No distinction of the flavours of produced quarks is made, and therefore the overall 
hadronic cross-section is measured, and is parameterised in terms of the hadronic width 
given by the sum over all quark final states, 

Thad ^ ; "^qq " 

The decays of the Z to neutrinos are invisible in the detectors and give rise to the 
"invisible width", r inv = N U T UU , where N v is the number of light neutrino species. The 
invisible width can be determined from the measurements of the decay widths to all visible 
final states and the total width, which is given by the sum over all partial widths, 

Tz = Tee + r w + r rr + Thad + Tinv • 

Because the measured cross-sections depend on products of the partial widths and also 
on the total width, the widths constitute a highly correlated parameter set. In order 
to reduce correlations among the fit parameters an experimentally-motivated set of six 
parameters is used to describe the total hadronic and leptonic cross-sections around the Z 
peak. These are 

• the mass of the Z, mz, and the total width, Tz] 
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• the "hadronic pole cross-section" , 

ee-Thad 
a h = —2 f2 ' 

m z i z 

• the ratios 

-R e = rhad/r ee , Rp = Thad/r w and R T = rh a d/r rT . 

The leading contribution from 7-Z interference is proportional to the product of the vector 
couplings of the initial and final states and vanishes at y/s = mz, but becomes noticeable 
at off-peak energies and therefore affects the Z mass. Because a determination of all quark 
couplings is not possible, the 7-Z interference term in the hadronic final state is fixed to 
its SM value. The implications of this are discussed in Section|Q|. 

Three additional parameters are needed to describe the leptonic forward-backward 
asymmetries for the processes e + e~ — > e + e~, e + e~ — ► /^ + /U~ and e + e~ — > t + t~. These 
are 

• the "pole asymmetries", A^b, an d -Ape- 

Contrary to the partial widths, the pole asymmetries are defined purely in terms of the 
real parts of the effective Z couplings, 

4{ = jA^k, with A,= — 

4 1 + 

Due to the smallness of the leptonic forward-backward asymmetry at y/s = mz, QED 
corrections are as large as .Ap B itself. The product of the axial vector couplings of the 
initial and final states determines the leading contribution of the 7-Z interference. This 
can be fixed with sufficient precision together with the vector couplings from simultaneous 
fits to the measured forward-backward asymmetries and cross-sections, requiring SM input 
only for the imaginary parts of the couplings. 

Differences between the pseudo-observables and the QED de-convoluted observables at 
y/s = mz, arising from the interference between photon and Z diagrams and from the 
interplay between the real and imaginary parts of the photon and Z couplings, are small 
in the SM. 0g eak , given in terms of the partial decay widths, agrees to better than 0.5 per- 
mit for both hadrons and leptons with the QED de-convoluted cross-sections without the 
photon exchange contribution at y/s = mz- This is only a small fraction of the Lep- 
combined experimental error. The difference between J^ B and the QED de-convoluted 
forward-backward asymmetry at the peak amounts to 0.0013, which is slightly larger than 
the LEP-combined error on J^ B . It is therefore important to treat the imaginary parts 
correctly, however, the measurements are not sensitive to variations of the imaginary parts 
within their SM expectation. 

The pseudo-observables introduced above cannot be considered as truly model-indepen- 
dent, because imaginary parts of the couplings as well as the 7-Z interference in the hadronic 



12 



final state need to be fixed to their SM values. This leads to small "Standard Model rem- 
nants" in any attempted "model-independent" definition of the pseudo-observables. More 
details about the treatment of imaginary parts and SM remnants in the theory programs 



TOPAZO UJ] and ZFITTER |L5|] are given in Reference [Tq. These computer codes make 
available the best current knowledge of QED and electroweak corrections within the min- 
imal Standard Model and also provide a (quasi) model-independent approach based on 
the parametrisation of the cross-sections and forward-backward asymmetries, as described 
above. 



4 Experimental results 

The experimental results presented here have been slightly modified from those published 
by the experiments [0-0] in order to facilitate the combination procedure. The four dedi- 
cated sets of experimental results for the combination are summarised in Table |5|. 

All fits are based on versions 6.23 of ZFITTER and 4.4 of TOPAZO. The published 
Aleph results were derived using version 6.10 of ZFITTER, which did not yet contain 



the improved treatment of fermion pairs radiated from the initial state [[17]]. For the 
combination presented here, the ALEPH measurements were reanalysed using version 6.23 
of ZFITTER. 

Each experiment used the combined energy error matrix described above (Table|j). 
This makes a small difference at the level of 0.1 MeV on mi and Tz and their errors only 
for L3, where uncertainties arising from the modelling of the radio frequency cavities are 
largest. 

In the Bhabha final state, the s-t interference has a small dependence on the value of 
the Z mass. Although practically negligible for a single experiment, a consistent treatment 
becomes important for the combination. Despite some different choices in the publications 
of the individual analyses, all experiments evaluate the t, s-t channel correction at their 
own value of the Z mass for the results presented here. The resulting interdependencies 
between the Z mass and the parameters from the Bhabha final state are explicitely included 
in the error correlation coefficients between mz and R e or A^g. 

The Lep experiments agreed to use a standard set of parameters for the calculation of 
the Standard Model remnants in the theory programs. The important parameters are the 
Z mass, m z = 91.187 GeV, the Fermi constant, G F = (1.16637±0.00001) x 10~ 5 GeV~ 2 g 



the electromagnetic coupling constant, a' 5 '(mz) = 

1/128.877±O.O9O0 [|H, the strong cou 
pling constant, a s (m z ) = 0.119±0.002 |20| , the top quark mass, m t = 174.3±5.1 GeV [pi 



and finally the Higgs mass, mn, which was fixed to 150 GeV. The dependence of the fit re- 
sults arising from the uncertainties in these parameters is almost negligible, as is discussed 
in SectionB.4.3. 



3 ) a^(mz) is the electromagnetic coupling constant at the scale of the Z mass for five quark flavours 
value and err> 
0.02804 ±0.00065 



the value and error given correspond to a correction due to hadronic vacuum polarisation of Aaj^ = 
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correlations 
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-.002 
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^FB 
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.043 


.003 


-.002 
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xWdf = 158/166 _ ~ L3 
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-.343 
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Rn 
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R T 
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X 2 /N d{ = 155/194 Opal 
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20.811 ± 0.058 


.001 


.020 


.222 


.093 


1.00 








R T 
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0.0089 ± 0.0045 
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-.005 
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a0,M 
^FB 
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Table 5: Results on Z parameters and their correlation coefficients. 
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5 Common uncertainties 



Important common errors among the results from all Lep experiments arise from several 
sources. These are the calibration of the beam energy, the theoretical error on the cal- 
culation of the small-angle Bhabha cross-section used as the normalisation reaction for 
all cross-section measurements, the theoretical uncertainties in the t-channel and s-t in- 
terference contribution to the differential large-angle Bhabha cross-section, the theoretical 
uncertainties in the calculations of QED radiative effects and, finally, from small uncertain- 
ties in the parametrisation of the electroweak cross-section near the Z resonance in terms 
of the set of pseudo-observables the four collaborations agreed upon. These common errors 
are quantified in the following sub-sections and are used in the combination. 

Other sources of common errors may arise from the use of common Monte Carlo gener- 
ators and detector simulation programs. However, each group uses its own tuning proce- 
dures and event selections which best suit their detector, and therefore the related errors 
are treated as uncorrelated among the experiments. 



5.1 Common energy uncertainties 

For the purpose of combining the experimental results at the parameter level, the errors on 
the centre-of-mass energy of each individual cross-section or asymmetry measurement, as 
given in Table [|, need to be transformed into errors on the extracted pseudo-observables. 
The first step is to scale the energy errors by factors of 1 ± e, while maintaining the exper- 
imental errors fixed. Typical values of e used are between 5% and 20%. Performing the 
standard fits with these scaled errors generates two pseudo-observable covariance matrices, 
V±, from which the covariance matrix due to energy errors, Ve, can be separated from the 
other errors, V exp , using the relation (V±) = (1 ± e) 2 (Vs) + (V cxp ). The validity of this pro- 
cedure was verified using a data set restricted to the hadronic cross-section measurements 
of the years 1993-1995, which were combined both at the cross-section level and at the 
parameter level. 



Table |15| and [16| in the appendix show the energy errors on the pseudo-observables 
extracted from the individual experimental data sets. The estimated energy errors differ 
slightly depending on which experimental data set is used to derive them. Combinations 
may be attempted based on each of them, or on the average. The central values and 
errors of each of the averaged parameters agree well to within 5 % of the error on that 
average. It is therefore most appropriate to take the average of the error estimates over 
the experiments as the common energy errors, which are shown in Table 0. 



5.2 Common t-channel uncertainties 

The t-channel and s-t interference contributions are calculated in the Standard Model 
using the program ALIBABA piZ| . The s-t interference contribution to the t-channel 
correction in Bhabha final states depends on the value of the Z mass. For the purpose of 
this combination, all experiments parametrise the t and s-t contributions as a function of 
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-0.0003 
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-0.0003 
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0.0003 



Table 6: Common energy errors for nine-parameter fits. Values are given as the signed square 
root of the covariance matrix elements in the same units as in Table^; elements above the di- 
agonal have been omitted for simplicity. The anti- correlation between electron and muon or tau 
asymmetries arises from the different energy dependence of the electron asymmetry due to the 
t-channel contribution. 



mz. This allows the t, s-t correction to follow the determination of mz in the fits, which 
results in a correlation between mz and R e or Apg. The change in correlation coefficients 
introduced by explicitly taking into account the m z dependence of the t channel in the fits 
is about +10 % for the mi-Re correlation and —10 % for the mz-ApB one. These correlation 
coefficients take the changes in R e and A^g properly into account when mz takes its average 
value in the combination of the four experiments. 

The theoretical uncertainty on the t-channel correction is discussed in detail in Refer- 
ence^. The size of the uncertainty is typically 1.1 pb for the forward cross-section and 
0.3 pb for the backward cross-section and depends slightly on the acceptance cuts p4 



All collaborations incorporate the theory uncertainty as an additional error on the electron 
pair cross-section and asymmetry. In order to evaluate the common error from this source, 
each collaboration performed two fits, with and without the theory error, and the quadratic 
differences of the covariance matrix elements for R e and Apg are taken as an estimate of 
the common error. The unknown error correlation between energy points below and above 
the peak is included in the error estimates by adding in quadrature the observed shifts 
in mean values of R e and Apg when varying this correlation between —1 and +1. The 
t, s-t related errors estimated by individual experiments are summarised in Table |TT| in the 
appendix. Since these are all very similar, the average shown in Table ^ is taken as the 
common error matrix. 





Re 


A 0,e 


Re 


0.024 




A 0,e 
^TB 


-0.0054 


0.0014 



Table 7: Common errors arising from the t-channel and s-t interference contributions to the 
e + e~ final states, given as the signed square root of the covariance matrix elements. 
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5.3 Common luminosity uncertainties 



The four collaborations use similar techniques to measure the luminosity of their data sam- 
ples by counting the number of small-angle Bhabha scattered electrons. The experimental 
details of the four measurements differ sufficiently that no correlation is believed to exist 
in the experimental component of the luminosity errors. All four collaborations, however, 
use BHLUMI 4.04 |[25|| , the best available Monte Carlo generator for small-angle Bhabha 
scattering, to calculate the acceptance of their luminosity counters. Therefore significant 
correlations exist in the errors assigned to the scale of the measured cross-sections due to 
the uncertainty in this common theoretical calculation. 

This uncertainty is estimated to be 0.061 % without applying a correction for 



the production of light fermion pairs, which is not calculated in BHLUMI, and enters as 
a contribution to the estimated error. A recent calculation of the contribution of light 
pairs J27| has allowed Opal to explicitly correct for light pairs and reduce its theoretical 
luminosity uncertainty to 0.054%. This is taken as common with the errors of the other 
three experiments, who between themselves share a mutual common error of 0.061 %. 

These errors almost exclusively affect the hadronic pole cross-section, and contribute 
about half its total error after combination. The common luminosity error also introduces a 
small contribution to the covariance matrix element between Tz and a®. This was neglected 
in the common error tables given above, as it had no noticeable effect on the combined 
result. 



5.4 Common theory uncertainties 

An additional class of common theoretical errors arises from the approximations and spe- 
cial choices made in the fitting codes. These comprise contributions from QED radiative 
corrections, including initial-state pair radiation, and the parametrisation of the differen- 
tial cross-section around the Z resonance in terms of pseudo-observables defined precisely 
at the peak and for pure Z exchange only. In order to estimate the uncertainties from the 
parametrisation of the electroweak cross-sections near the Z resonance the two most ad- 
vanced calculational tools, ZFITTER [15] and TOPAZ0 [Hj, were compared. In addition, 



there are "parametric uncertainties" arising from parameters of the SM that are needed to 
fix the SM remnants. 



5.4.1 QED uncertainties 

The effects of initial state photon and fermion pair radiation lead to the large corrections in 
the vicinity of the Z resonance illustrated in Figure [l], and therefore play a central role in the 
extraction of the pseudo-observables from the measured cross-sections and asymmetries. 
Such large radiative effects have to be seen in contrast to the experimental precision, which 
is below the per-mil level in the case of the hadronic cross-section. 

The most up-to-date evaluations of photonic corrections include the leading contribu- 
tions up to 0(a 3 ). Two different schemes are available to estimate the remaining uncer- 
tainties: 
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1. KF: 0(...a 2 L 2 , a 2 L, a 2 L°) calculations [E8| including the exponentiation scheme of 
Kuraev-Fadin |j| with 0(a 3 L 3 ) gg. 



2. YFS: the 2 nd order inclusive exponentiation scheme of Reference |3lj [32], based on the 



YFS approach Third order terms are known and have only a small effect j34] 



Differences between these schemes, which are both implemented in ZFITTER, TOPAZO 



and MIZA [35], and uncertainties due to missing higher order corrections [34], amount to 



at most ±0.1 MeV on m z and T z , and ±0.01 % on cx°. 

The influence of the interference between initial and final state radiation on the ex- 
tracted parameters has also been studied recently []36| , and uncertainties on m z of at most 



±0.1 MeV from this source are expected for the experimental results given with only a small 
cut on s', the effective squared centre-of-mass energy after photon radiation from the initial 
state. The uncertainties due to the extrapolation of the leptonic s-channel cross-sections 
to full angular acceptance and from large to small s' are different among the experiments 
and are believed to be largely uncorrelated. 

QED related uncertainties are dominated by the radiation of fermion pairs from the 
initial state. Starting from the full second order pair radiator |28|, B7], a simultaneous 



exponentiation scheme for radiated photons and pairs was proposed in Reference BR A 



third-order pair radiator was calculated recently [17] and compared with the other ex- 
isting schemes, which are all available in the latest version of ZFITTER. Independent 
implementations of some schemes exist in TOPAZO and in MIZA. The largest uncertainty 
arises from the contribution of hadronic pairs. The maximum differences are 0.3 MeV on 
m z , 0.2 MeV on T z and 0.015 % on < 

In summary, comparing the different options for photonic and fermion pair radiation 
leads to error estimates of ±0.3 MeV on m z and ±0.2 MeV on T z . The observed differences 
in a® are slightly smaller than the error estimate of ±0.02 % in Reference^, which is 
therefore taken as the error for QED uncertainties. 

5.4.2 Parametrisation of line shape and asymmetries 



In a recent detailed comparison of TOPAZO and ZFITTER ||16|| , cross-sections from Stan- 
dard Model calculations and from the model-independent parametrisation were considered. 
Uncertainties on the pseudo-observables arise from differing choices in the parametrisa- 
tion of the electroweak cross-sections near the Z resonance. In order to determine these 
TOPAZ0-ZFITTER differences, each of the two codes have been used. For practical 
reasons, cross-sections and forward-backward asymmetries were calculated with TOPAZO 
and then fitted with ZFITTER. Errors were assigned to the calculated cross-sections and 
forward-backward asymmetries which reflected the integrated luminosity taken at each 
energy, thus ensuring that each energy point entered with the appropriate weight. 

The dominant part of the small differences between the two codes results from details of 
the implementation of the cross-section parametrisation in terms of the pseudo-observables. 



4 ) Third-order terms for the KF scheme had also been calculated earlier [filf 
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0.0001 


0.0001 


0.001 


0.004 


0.0001 



Table 8: Differences in fit results obtained with TOPAZ0 and ZFITTER, taken as common 
systematic errors. 

This is particularly visible for the off-peak points, where the assignment of higher-order 
corrections to the Z resonance or to the SM remnants is not in all cases unambiguous. 
The size of the differences also depends on the particular values of the pseudo-observables, 
since these do not neccessarily respect the exact SM relations. Slightly different choices 
are made in the two codes if the SM relations between the pseudo-observables are not 
fulfilled. Finally, variations of factorisation schemes and other options in the electroweak 
calculations may affect the fit results through the SM remnants, but were found to have a 
negligible effect. 

In Table |] differences between TOPAZ0 and ZFITTER are shown, which are taken 
as common systematic errors. They were evaluated around the set of pseudo-observables 
representing the average of the four experiments, where the cross-sections and asymmetries 
were calculated for full acceptance with only a cut on s 1 > 0.01 s. The only important 
systematic error of this kind is the one on Re, which amounts to 15 % of the combined 
error. 

Putting all sources together, overall theoretical errors as listed in Table || are obtained, 
and these are used in the combination. 
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Table 9: Theoretical errors from fit programs, i. e. photon and fermion-pair radiation and model- 
independent parametrisation, given as the signed square root of the covariance matrix elements. 
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5.4.3 Parametric uncertainties 

Through the SM remnants the fit results depend slightly on the values of some SM pa- 
rameters. Varying these within their present experimental errors, or between 100 GeV and 
1000 GeV in case of the Higgs boson mass, leads to observable effects only on the Z mass, 
which is affected through the 7-Z interference term. The dominant dependence is on mH, 
followed by a' 5 '(mz). 

The effect on mz from a variation of l/a' 5 '(mz) by its error of ±0.090 is +0.05 MeV, 
which is negligibly small compared to the systematic error on mz arising from other QED- 
related uncertainties (see Tabled). The change in mz due to m H amounts to +0.23 MeV 
per unit change in log 10 (mH/GeV). Note that this is small compared to the total error on 
mz of ±2.1 MeV and is not considered as an error, but rather as a correction to be applied 
once the Higgs boson mass is known. The consequences of a completely model-independent 
treatment of the 7-Z interference in the hadronic channel are discussed in Section 



6 Combination of results 

The combination of results on the Z parameters is based on the four sets of results on the 
nine parameters m z , T z , cr£j, R e , R T , Apg, Apg and Apg and the common errors given 
in the previous chapter. 

For this purpose it is necessary to construct the full (4x9) x (4 x 9) covariance matrix 
of the errors. The four diagonal 9x9 matrices consist of the four error matrices specified 
by each experiment (Table|J). The 9x9 common error matrices build the off-diagonal 
elements. Some theoretical uncertainties must also be added to the diagonal matrices, 
since they are not contained in the individual experimental matrices. 

A symbolic representation of this matrix is shown in Table [10|. Each table element 
represents a 9 x 9 matrix; (C exp ) for exp = A, D, L and O are the covariance matrices 
of the experiments (Tabled), and (C c ) = {Ce) + {Cc) + {Ct) + (Cqeda) is the matrix of 
common errors. Here {Ce) (Table ^) is the error matrix due to Lep energy uncertainties, 
{Cc) (Section |5.3D arises from the theoretical error on the small-angle Bhabha cross-section 
calculations, {Ct) (Table|7|) contains the errors from the t-channel treatment in the e + e~ 
final state, and (CQED.th) contains the errors from initial state photon and fermion pair 
radiation and from ambiguities in the model- independent parametrisation (Table ^). Since 
the latter errors were not included in the experimental error matrices, they were also added 
to the block matrices in the diagonal of Table [TO]. 

The combined parameter set and its covariance matrix are obtained from a \ 2 minimi- 
sation, with 

x a = (x-x m ) T (c)- 1 (x-x 1II ), 

where (X — X m ) is the vector of residuals of the combined parameter set to the individual 
results. 

Some checks of the combination procedure outlined above are described in the following 



subsections, and the combined results are given in the tables of Section IO. 
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Table 10: Symbolic representation of the covariance matrix, (C), used to combine the line 
shape and asymmetry results of the four experiments. Elements above the diagonal are the 
same as those below and are left blank for simplicity. The components of the matrix are 
explained in the text. 



6.1 Multiple-mz fits 

In 1993 and 1995, the two years when Lep performed precision scans to measure the Z line 
shape, the experimental errors are very comparable, but the Lep energy was appreciably 
better understood in 1995 than in 1993. In determining the optimum value of mz, therefore, 
the four experiments combined should give more weight to the 1995 data than they each 
do in their independent determinations. To quantify this issue the measurements of each 
experiment were fitted to determine independent values of mz for the periods 1990-1992, 
1993-94 and 1995. In this "eleven-parameter fit", each of the three mass values m| 0-92 , 
m| 3-94 and m| 5 has its specific energy error reflecting the different systematic errors on the 
absolute energy scale of Lep. In the combination, the relative importance of energy related 
and independent experimental errors on the mass values is properly treated. The input 
and the common energy errors, estimated in the same way as for the nine parameters, are 
shown in Section A- 3. 



When the three values of mz are condensed into a single one, the effects of the time 
dependence of the precision in the energy calibration is taken into account. The difference 
of —0.2 MeV w.r.t. the mz value from the nine-parameter fits corresponds to 10% of the 
combined error. All other parameters are identical to their values from the nine-parameter 
fit to within less than 5 % of the combined error. This result justifies using the standard 
combination based on the nine parameters. Tables |T5| and [16] in the appendix show the 
energy errors on the pseudo-observables estimated from each individual experimental data 
set. The estimated energy errors differ slightly depending on which experimental data set 
is used to derive them. Combinations may be attempted based on each of them, or on the 
average. The central values and errors of each of the averaged parameters agree well within 
5 % of the error on that average. It is therefore most appropriate to take the average of 
the error estimates over the experiments as the common energy errors, which are shown in 
Table |. 

The averages over the four experiments of the three values m^ 0-92 , m| 3-94 and m| 5 also 
provide a cross-check on the consistency of the energy calibration, which dominates the 
errors on mz in each of the periods considered. Using the energy errors of Table |T9| allows 
the correlated and uncorrelated parts of the errors on the mass differences to be quantified. 
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This is shown in Figure |3|. The differences between these values for the Z mass amount 
to K°- 92 - ml 3 " 94 1 = 31%, K " 92 - mf\ = 56% and |m 93 " 94 - mf\ = 43% of the 
uncorrelated error, i. e. the three Z mass values are well consistent. 



1990-1992 

91.1904±0.0065 
1993-1994 

91.1882±0.0033 
1995 

91.1866±0.0024 
average 

91.1874±0.0021 



91.185 91.19 91.195 
m z [GeV] 



Figure 3: mz in GeV for different periods of data taking, before 1993, 1993-1994 and 
1995. The second, smaller error bar represents the correlated error component of 1.2MeV 
between m| 3-94 and m| 5 . m 9 / 1-92 is essentially uncorrelated with the other two. 



6.2 Fits with lepton universality 

All experiments provide fits to their measured cross- sect ions and asymmetries with lepton 
universality imposed, i. e. R e , and R T are replaced by Rg and A^g, A^'g and Ap^ 
get replaced by A^' B in the model-independent parametrisation of the differential cross- 
section. Here Rn is defined for massless leptons. The individual experimental results and 
the correlation matrices are given in Table|Tg| in the appendix. 

Comparing these five-parameter results with the nine-parameter results of Table [|, there 
is a noticeable change in mz of a few tenth of MeV in all experiments. This is a consequence 



of the dependence of the t-channel correction on mz, as discussed in Section [5T2| . When 
R e and Ap^ are replaced by the leptonic quantities Rt and A^g, their correlation with the 
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Z mass leads to a shift, which is driven by the (statistical) difference between R e and Ri 
and and Ap B . Similarly, replacing i? e and Apg from the values of a single experiment 
by the Lep average introduces a shift in mz in the presence of these particular correlation 
coefficients. Such shifts should become smaller when averaged over the four experiments. 
Indeed, the average of the shifts is only — 0.2MeV. 

Another aspect of the five-parameter results concerns the role of the common t-channel 
uncertainty in the averages over the leptonic measurements to obtain Rg and Apg. If the 
average over the leptonic measurements is performed by each experiment individually, the 
weight given to the electron channel is larger than for the case where the averages over 
individual lepton species are averaged at the end. In the latter procedure the weight of the 
electrons relative to the muon and tau final states is reduced due to the common t-channel 
error. Extracting the results with lepton universality from the nine parameters is therefore 
the appropriate method. 



6.3 Shifts for halved experimental errors 

When the average over the experiments is performed at the parameter level, information 
on the individual contribution of particular data points to the average is lost. Performing 
the average over the data points instead may therefore lead to changes of the relative 
importance of independent experimental errors w.r.t. the common errors. The examples 
of mz and the importance of the t-channel errors for Ri, as discussed in the previous 
subsections, provide good illustrations of such effects. It was demonstrated that averaging 
the shifts in mz which each experiment observed when halving its experimental errors to 
simulate the generic "combined" experiment also reproduced the results of the full fit to 
the combined hadronic cross-section measurements. 

While mz is properly treated by the eleven-parameter fits, other pseudo-observables 
may suffer from similar changes due to shifts of the weights. Changes in central values 
when halving the independent experimental errors in each experiment can be used as a 
monitoring tool for these parameters as well. The average of these changes over the four 
Lep experiments serves to control the differences between an average at the parameter level 
compared to a full cross-section average. Of course, this assumes that all measurements 
from individual experiments enter into the average with the same weight. The observed 
shifts are summarised in Table 0. The shift downwards in mz of 0.3 MeV corresponds to 
the slightly smaller shift of 0.2 MeV already seen in the multiple-mz fits. 

Thus, the average changes in mz, 0"°, -R e , A^'g and A^'g amount to about 10% of the 
combined errors, in all other cases they are even smaller. 

This cross-check provides an estimate of the magnitude of the changes in the final 
results that would arise from a combination at the cross-section level. The averaging at 
the parameter level is equal to this within at most 10% of the combined error. 
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A 


D 


L 


o 


average 


% of error 


m z [GeV] 


-0.0006 


0.0000 


-0.0004 


-0.0001 


-0.00028 


13 


T z [GeV] 


-0.0002 


+0.0001 


-0.0004 


0.0000 


-0.00013 


5 


< [nb] 


+0.006 


0.000 


+0.008 


+0.0036 


+0.0037 


10 


R e 


+0.004 


+0.017 


0.000 


+0.004 


+0.0063 


13 


Rfj, 


0.000 


0.000 


0.000 


+0.001 


0.0000 





R r 


0.000 


0.000 


-0.001 


+0.002 


+0.0003 


1 


A 0,e 
■"FB 


-0.0001 


-0.0003 


0.0000 


-0.0000 


-0.00011 


5 


-"FB 


+0.0002 


+0.0003 


0.0000 


+0.0001 


+0.00014 


11 


A 0,r 

-^fb 


+0.0002 


+0.0003 


0.0000 


+0.0001 


+0.00015 


9 



Table 11: Shifts in central values of the fitted pseudo-observables seen when halving the indepen- 
dent experimental errors, for individual experiments and average. 

6.4 Influence of the 7-Z interference term 

In the nine-parameter analyses discussed here, the 7-Z interference terms in the differen- 
tial cross-sections for leptons are expressed using the effective coupling constants and the 
electric charges of the electron, Q e , and the final state fermion, Q f (see equation for the 
differential cross-section in Section^. For the hadronic final state, however, the 7-Z inter- 
ference terms are fixed to the SM values, as individual quark flavours are not separated. 
Fits with a free interference term are possible in the S-Matrix scheme ||39|| . The OPAL 
collaboration also studied a different approach |4[] based on an extension of the standard 
parameter set. In the S-Matrix approach the interference terms are considered as free and 
independent parameters. The hadronic interference term is described by the parameter 
Jtot d 5 given in the SM by 

V27ra(m z ) q 

Note that the running of a as well as final state QED and QCD corrections are also included 
in the definition of the S-Matrix parameters. The SM value of j t h ^ d is 0.21 ± 0.01. 

The dependence of the nine parameters on the hadronic 7-Z interference term is studied 
here by considering a set of ten parameters consisting of the standard nine parameters 
extended by the parameter j t h ^. d from the S-Matrix approach. The 7-Z interference terms 
in the lepton channels are fixed by the leptonic Z couplings. As already observed in S- 
Matrix analyses of the LepI data 0, ||, a large anti-correlation between m z and j t h ^. d 
appears, leading to errors on m z enlarged by a factor of almost three. The studies show 
that the dependence of m z on is given by 

= -1.6MeV/0.1. 

dj t h o a t d ' 

The changes in all other parameters are below 20 % of the errors in the combination for a 
change in j t h ^ d of 0.1 . 
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Better experimental constraints on the hadronic interference term are obtained by in- 
cluding measurements of the hadronic total cross-section at centre-of-mass energies further 
away from the Z pole than just the off-peak energies at Lep I. Including the measurements 
of the TRISTAN collaborations TOPAZ gg and VENUS [|TJ at ^/(s)=58GeV, the error 
on j^f is reduced to ±0.1, while its central value is in good agreement with the SM ex- 
pectation. Measurements at centre-of-mass energies above the Z resonance at Lep II also 
provide constraints on j t h ^. d , but in addition test modifications to the interference terms 



arising from the possible existence of a heavy Z' boson p2|-|45[ . 

The available experimental constraints on j t h ^ d thus lead to uncertainties on m z , inde- 
pendent of SM assumptions in the hadronic channel, which are already smaller than its 
error. No additional error is assigned to the standard nine-parameter results from effects 
which might arise from a non-SM behaviour of the 7-Z interference. 



6.5 Direct Standard Model fits to the measured cross-sections 
and forward-backwards asymmetries 

Since an important use of the combined results presented here is to determine parameters 
of the SM and to make tests on its validity, it is crucial to verify that the parameter set 
we have chosen for the combination represents the four sets of experimental measurements 
from which they were extracted in a manner adequate to this purpose. When the set 
of pseudo-observables is used in the framework of the Standard Model, the role of m z 
changes from an independent parameter to that of a Lagrangian parameter of the theory. 
This imposes additional constraints which can be expected to shift the value of m z . 

To check whether the nine parameters adequately describe the reaction to these con- 
straints, each collaboration provided results from direct Standard Model fits to their 
cross-section and asymmetry data. The comparison of these results with those obtained 



by using the set of pseudo-observables as fit input is shown in Table |12j. mn and a s 
were free parameters in these fits, while the additional inputs m t =174.3±5.1 GeV and 
Aa£2i = 0.02804 ± 0.00065 (corresponding to l/a^(m z ) = 128.877 ± 0.090) provided 
external constraints. 

Most noticeably, significant shifts in m z are observed in some individual data sets, which 
cancel to almost zero in the average over the four experiments. One anticipated source 
of these shifts was already mentioned: the Z couplings defining the 7-Z interference term 
depend on toh, which is allowed to move freely in the fit for SM parameters, but is fixed 
to 150 GeV in the model-independent fit for the extraction of the pseudo-observables. The 
approximate values of mn preferred by the SM fit to the cross-sections and asymmetries 
are indicated in the second part of the table. Using the dependence of m z on the value of 



mH given in Section |5.4.3| , the differences in m z can be corrected to a common value of the 
Higgs mass of mH=150 GeV, as is shown in the last line of Table fL2[ Hence the net average 
difference in m z directly from the data or through the intermediary of the nine parameters 
is less than 0.1 MeV. Shifts in the other SM parameters, in the individual data sets as well 
as in the average, are all well under 5 % of the errors, and therefore also negligible. 
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A 


D 


L 





average 


% of error 


xV^Vdf 


174/180 


184/172 


168/170 


161/198 




Am z [MeV] 


-0.7 


+0.5 


0.0 


+0.1 


-0.03 


1 


Am t [GeV] 


0.0 


0.0 


0.0 


0.0 


0.0 


<2 


Alog 10 (m H /GeV) 


-0.01 


+0.04 


+0.02 


+0.04 


+0.02 


4 


Aa s 


0.0000 


-0.0002 


+0.0002 


+0.0002 


+0.0001 


4 




+0.00002 


-0.00004 


0.00000 


-0.00004 


-0.00002 


2 


fit value 














of m H [GeV] 


40. 


10. 


35. 


390. 






Am z [MeV] 














corrected to 














m H =150GeV 


-0.6 


+0.7 


+0.1 


0.0 


+0.05 


2 



Table 12: Shifts in SM parameters, from direct SM fit to the cross-sections and forward-backward 
asymmetries w.r.t. fits to the nine-parameter results. The numbers in the lowest part of the table 
give the shifts in m z if the results from the first line are corrected to a common value of the Higgs 
mass of 150 GeV. 



The conclusion of this study is that SM parameters extracted from the pseudo-ob- 
servables are almost identical to the ones that would be extracted from the combined 
cross-sections and asymmetries. Within the SM the combined set of pseudo-observables 
provides a description of the measurements of the Z parameters that is equivalent to the 
full set of cross-sections and asymmetries. This is also true for any theory beyond the SM 
which leads to corrections that are absorbed in the pseudo-observables. An exception to 
this are those theories with an additional Z'-bosons which have significant modifications of 
the 7-Z interference term. (See the discussion in Section |£T4} ) 



6.6 Combined results 

The result of the combination of the four sets of nine pseudo-observables of Tablej^, includ- 



ing the experimental and common error matrices shown in Table [10|, is given in Table |i3[ 
The value of the Higgs boson mass was assumed to be 150 GeV and is relevant only for 
the value of mz, which changes by +0.23 MeV per unit change in log 10 (mH/GeV). (See 
Section jXj .) 

The value of y 2 per degree of freedom of the combination is 32.6/27 and corresponds 
to a probability of 21 % to find an agreement among the four sets of measurements which 
is worse than the one actually observed. The correlation matrix of the combined result 
shows significant correlations of cr° with T z and R e , i? M and R T and between R e and A^^. 

A comparison of the leptonic quantites R e , i? M and R T and of Ap£, Apg and Apg shows 
that they agree within errors. Note that R T is expected to be larger by 0.23% because of 
r mass effects. Figure |] shows the corresponding 68% level contours in the Rf-J^ B plane. 

Imposing the additional requirement of lepton universality in the combination leads to 
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without lepton universality 


correlations 


xV^Vdf = 32.6/27 


mz 


r z 


< 


R c 




R T 


a U, e a u, fj, 
^FB ^FB 


A u ' 
^TB 


m z [GeV] 


91.1876± 0.0021 


1.00 
















r z [GeV] 


2.4952 ± 0.0023 


-.024 


1.00 














< [nb] 


41.541 ± 0.037 


-.044 


-.297 


1.00 












R c 


20.804 ± 0.050 


.078 


-.011 


.105 


1.00 










Rfi 


20.785 ± 0.033 


.000 


.008 


.131 


.069 


1.00 








R T 


20.764 ± 0.045 


.002 


.006 


.092 


.046 


.069 


1.00 






A°' e 
•"FB 


0.0145 ± 0.0025 


-.014 


.007 


.001 


-.371 


.001 


.003 


1.00 




0.0169 ± 0.0013 


.046 


.002 


.003 


.020 


.012 


.001 


-.024 1.00 




A 0,r 
^FB 


0.0188 ± 0.0017 


.035 


.001 


.002 


.013 


-.003 


.009 


-.020 .046 


1.00 



with lepton universality 



X 2 /^df = 36.5/31 


mz 


r z 


< 


Re 


A^ 


m z [GeV] 


91.1875± 0.0021 


1.00 










T z [GeV] 


2.4952 ± 0.0023 


-.023 


1.00 








< [nb] 


41.540 ± 0.037 


-.045 


-.297 


1.00 






R £ 


20.767 ± 0.025 


.033 


.004 


.183 


1.00 




*TB 


0.0171 ± 0.0010 


.055 


.003 


.006 


-.056 


1.00 



Table 13: Result of the combination of the four sets of nine pseudo-observables from Table^. 



the results shown in the second part of Table 13. Here Re is not a simple average over 



the three lepton species, but refers to Z decays into pairs of one massless charged lepton 
species. The value of x 2 /N^ of 36.5/31 for the combination of the four sets of nine pseudo- 
observables into the five parameters of Table |13"| corresponds to a x 2 -probability of 23%. 
The central ellipse in Figure |] shows the 68 %-CL contour for the combined values of Ri 
and Ap B determined from all three lepton species. 

In principle, the average over the four experiments can also be performed at the level of 
the five parameter results of Section |BT2] . When this is attempted, good agreement is seen 
with the results in the last line of Table |T3|, except for R £ , where the difference amounts to 
0.005 or 20 % of the total error. The origin of this shift is the common t-channel error, as 
discussed in Section|673 above. 



6.7 Parameter Transformations 

Additional pseudo-observables, more familiar than the experimentally motivated set of 
Table |i~3"l, can be obtained by parameter transformations. The partial Z decay widths are 
summarised in Table [T3|; it should be noted that these have larger correlations than the 
original set of results. If lepton universality is imposed, the value of r had also changes and 
its error is reduced, because T ee in the relation between the hadronic pole cross-section 
and the partial widths is replaced by the more precise value of T«. The invisible width, 
Tinv = r z — Thad - r ee — r w - r rr , is also shown in the table. The leptonic pole cross-section, 
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0.022 



0.018 
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0.01 
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Figure 4: Contour lines (68 % CL) in the Rg -A^b plane for e + e~, and t + t~ final states 

and for all leptons combined. For better comparison the results for the r lepton are corrected 
to correspond to the massless case. The SM prediction for mz=91.1875G(N , m t =i7^.5GeV ; 
mH=500GeV and a s =0.119 is also shown. The lines with arrows correspond to the variation 
of the SM prediction when nit, m n and a s are varied in the intervals mt = 174.3 ± 5.1 GeV, 



m H 



300 



+700 
-200 



GeV, and a s = 0.119 ± 0.002, respectively. The arrows point in the direction of 



increasing values ofm t , ran and a s 



a° tl defined as 



12vr Tl 



Or 



rrin 



F 2 
1 z 



in analogy to cr{J, is given in the last line of Table [H|. Because QCD final state corrections 
appear quadratically in the denominator via Tz, a° t has a higher sensitivity to a s than 
or Re, where the dependence on QCD corrections is only linear. 

Assuming only standard particles, the invisible width is compatible with the SM hy- 
pothesis of decays into the three neutrino species, r inv = 3T UU . The "number of neutrinos" , 
N u , is calculated according to 



U / SM 



The Standard Model value for the ratio of the partial widths to neutrinos and to charged 
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without lepton universality 


correlations 


r had [MeV] 1745.8 ±2.7 
r ee [MeV] 83.92±0.12 
T w [MeV] 83.99±0.18 
T TT [MeV] 84.08±0.22 


1.00 
-0.29 1.00 
0.66-0.20 1.00 
0.54-0.17 0.39 1.00 


with lepton universality 


T inv [MeV] 499.0 ±1.5 
r had [MeV] 1744.4 ±2.0 
T u [MeV] 83.984±0.086 


1.00 
-0.29 1.00 
0.49 0.39 1.00 


r inv /r« 5.942±0.016 




a° e [nb] 2.0003±0.0027 





Table 14: Partial Z decay widths and correlation coefficients. 

leptons is 1.9912±0.0012, where the uncertainty arises from variations of the top quark mass 
within its experimental error and of the Higgs mass within 100 GeV < mn < 1000 GeV. 
With the measured value of r inv / T u = 5.942 ± 0.016, the number of light neutrino species 
is determined to be 

N v = 2.9841 ± 0.0083 . 

This may also be turned into a quantitative limit on extra, non-standard contributions to 
the invisible width, i. e. not originating from Z — > vT>, by taking the difference between the 
value given in Table |i~4| and the Standard Model expectation of (r inv ) SM = 501.7ta9 MeV. 
This gives Ar inv = -2.7±\'l MeV, or expressed as a limit, Ar inv < 2.0 MeV @ 95% CL; 
here, the limit was conservatively calculated allowing only positive values of r^ v . 

The effective axial- vector and vector couplings of the Z to leptons, g l A and g{r, may be 
expressed in terms of the effective Veltman p parameter Kl and the effective weak mixing 
angle, sin 2 ^, by 



A 



9v 



3 • 



pT(4 



where J| = — ^ is the weak iso-spin of charged leptons. 



The leptonic pole asymmetry, Ap B , depends only on the ratio of the universal lepton 
couplings Qy I g e A , as is easily seen from its definition given in Section^], and thus directly 
determines the effective weak mixing angle, 



sin 



2 /ilept 
^cff 




if the sign of the ratio of couplings is chosen to be positive, in agreement with studies of 
the polarisation of r leptons at Lep. 
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The leptonic width is almost entirely given by the parameter p 1 ^, with a small contri- 
bution from the weak mixing angle entering through the vector couplings. A simultaneous 
fit of Tee and Ap B for p 1 ^ and sin 2 ^^ results in 

p|f ff pt = 1.0048 ±0.0011, 
sin 2 #' c ff pt = 0.23099 ± 0.00053 , 

with an error correlation coefficient of 27% . 

Information on the effective p parameter contained in r; nv and I^ad does not significantly 
improve on this result. From r; nv , the effective p parameter for neutrinos, p v eS = (2g u ) 2 , 
is determined to be Pg ff = 1.0027 ± 0.0030. The extraction of p from the hadronic width 
would require external constraints on the strong coupling constant and SM assumptions 
on flavour-dependent corrections. 

The result on sin 2 ^^* derived here from the leptonic forward-backward asymmetries 
agrees well with the value recently published by SLD based on measurements of the left- 



right polarisation asymmetry of hadron production at the Z resonance [47]. The above 
value of pjf* is 4.4 standard deviations greater than its tree-level value of one, in good 
agreement with full SM calculations, and thus clearly shows the presence of genuine elec- 
troweak corrections. 



7 Conclusion 

The combination procedure adopted by the LEP electroweak working group for the four 
sets of pseudo-observables derived by the Lep experiments from the measured hadronic and 
leptonic cross-sections and the leptonic forward-backward asymmetries at centre-of-mass 
energies around the Z resonance has been described. 

The combination procedure averages parameters extracted from the measurements of 
each individual experiment. This approximates a statistically optimal average that would 
be based on the measured cross-sections and forward-backward asymmetries to better than 
10 % of the combined errors, i. e. the chosen set of parameters adequately represents the 
full set of measurements. The technical precision of the adopted combination procedure is 
around 5 % of the combined errors. Using each of the four sets of input parameters in the 
framework of the minimal Standard Model yields results which are on average identical, 
to within 5 % of the combined errors, to those obtained directly from Standard-Model fits 
to the measured cross-sections and asymmetries of each individual experiment. 

Detailed studies of common systematic errors were performed. The dominant contri- 
bution of ±1.7 MeV to the combined error on the Z mass arises from the calibration of the 
energy of the beams in Lep. The dominant contribution of ±0.025 nb to the uncertainty on 
the hadronic cross-section at the pole of the Z resonance arises from the theoretical error 
on the small-angle Bhabha cross-section. The errors on all other parameters are dominated 
by independent experimental or statistical errors. 
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The combined Lep results on the mass and width of the Z, on the hadronic pole 
cross-section, on the ratio of the hadronic and leptonic partial width and on the pole 
forward-backward asymmetry are 

m z = 91.1875 ± 0.0021 GeV, 

T z = 2.4952 ±0.0023 GeV, 

al = 41.540 ± 0.037 nb, 

Re = 20.767 ±0.025, 

A^'b = 0.0171 ±0.0010. 

Correlations are typically small (< 5%), being significant only between a{J and Tz 
(—30%) and between 0"° and Re (18%). The full set of results and the error correlation 



matrices are shown in Table 113. 
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A Appendix 



A.l Common errors estimated by individual experiments 







mz 


Tz 


< 


Re 


Aleph 


m z [GeV] 
T z [GeV] 


0.0017 
-0.0004 


0.0012 








< [nb] 


-0.0028 


-0.0024 


0.011 






R c 


0.0013 


-0.0015 


0.007 


0.012 


Delphi 


m z [GeV] 
T z [GeV] 


0.0016 
-0.0005 


0.0012 








< [nb] 


-0.0025 


-0.0024 


0.009 






Re 


0.0014 


0.0000 


0.004 


0.016 


L3 


m z [GeV] 
Tz [GeV] 


0.0016 
-0.0008 


0.0013 








< [nb] 


-0.0009 


-0.0030 


0.011 






R c 


0.0020 


-0.0021 


0.010 


0.013 


Opal 


m z [GeV] 
T z [GeV] 


0.0017 
-0.0005 


0.0012 








[nb] 


-0.0011 


-0.0028 


0.011 






Re 


0.0019 


-0.0019 


0.008 


0.012 



Table 15: Common energy errors, from individual data sets for the nine-parameter fits. Values 
are given as the signed square root of the covariance matrix elements; elements above the diagonal 
have been omitted for simplicity. 
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A 0,c 


A 0,M 


A°' T 


Aleph 


A U,e 


0.0003 








/\0,M 


-0.0003 


0.0003 






A 0,r 
^B 


-0.0003 


0.0003 


0.0003 


Delphi 


A U,c 


0.0004 








/\0,M 
^B 


-0.0003 


0.0003 






A 0,r 
^B 


-0.0003 


0.0003 


0.0003 


L3 


A U,c 
^B 


0.0003 








^B 


-0.0002 


0.0003 






^B 


-0.0002 


0.0002 


0.0003 


Opal 


A U,c 
"^FB 


0.0004 








^B 


-0.0003 


0.0003 






^B 


-0.0003 


0.0003 


0.0003 



Table 16: Common energy errors for forward-backward asymmetries, estimated from the data 
sets of individual experiments. 







Re 


^B 


Aleph 


R e 


0.025 






A°' e 


-0.0056 


0.0013 


Delphi 


R e 


0.025 






A°' c 

•™FB 


-0.0058 


0.0016 


L3 


Re 


0.021 






A°' c 

™FB 


-0.0046 


0.0010 


Opal 


R c 


0.025 






^B 


-0.0057 


0.0015 



Table 17: Estimates of the t-channel errors of individual experiments. 
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A. 2 Results with lepton universality 



correlations 



X 2 /N di = 172/180 Aleph 



m z [GeV] 91.1893 ± 0.0031 
T z [GeV] 2.4959 ± 0.0043 
erg [nb] 41.559 ± 0.057 
R e 20.729 ± 0.039 
A^'b 0.0173 ± 0.0016 


1.00 

.038 1.00 
-.092 -.383 1.00 
.033 .011 .246 1.00 
.071 .002 .001 -.076 1.00 


x 2 /AT df = 183/172 Delphi 


m z [GeV] 91.1863 ± 0.0028 
T z [GeV] 2.4876 ± 0.0041 
erg [nb] 41.578 ± 0.069 
R e 20.730 ± 0.060 
A^'b 0.0187 ± 0.0019 


1.00 

.046 1.00 
-.070 -.270 1.00 
.028 -.006 .242 1.00 
.095 .006 -.005 .0001.00 


X 2 /N di = 163/170 L3 


m z [GeV] 91.1894 ± 0.0030 
T z [GeV] 2.5025 ± 0.0041 
erg [nb] 41.536 ± 0.055 
20.809 ± 0.060 
A^'b 0.0192 ± 0.0024 


1.00 

.068 1.00 

.014 -.348 1.00 

.067 .020 .111 1.00 

.041 .020 .005 -.024 1.00 


X 2 /N di = 158/198 Opal 


m z [GeV] 91.1853 ± 0.0029 
T z [GeV] 2.4947 ± 0.0041 
erg [nb] 41.502 ± 0.055 
Re 20.822 ± 0.044 
A^-b 0.0145 ± 0.0017 


1.00 

.051 1.00 

.030 -.352 1.00 

.043 .024 .290 1.00 

.075 -.005 .013 -.017 1.00 



Table 18: Results on Z parameters and error correlation matrices by the four experiments with 
lepton universality imposed. 
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A. 3 Results from eleven-parameter fits 



Aleph 



r z 


[GeV] 


2 


4957 


± 





0043 


< [nb] 


41 


559 


± 





057 




Re 


20 


694 


± 





075 




Rn 


20 


801 


± 





056 




R T 


20 


709 


± 





062 




»0,e 





0184 


± 





0034 




"^FB 





0173 


± 





0025 




A 0,r 





0171 


± 





0028 




[GeV] 


91 


1928 


± 





0092 


mi 3 " 94 


[GeV] 


91 


1926 


± 





0046 


95 


[GeV] 


91 


1852 


± 





0043 



correlation matrix 



1.000 


-0.384 





001 





012 





003 





001 





000 





000 





015 


-0 


Oil 





055 


-0.384 


1.000 





137 





167 





152 


-0 


004 





003 





003 


-0 


043 


-0 


037 


-0 
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-0 
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-0 
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1 
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Delphi 
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[GeV] 


2.4878 


± 





0041 




< [nb] 


41.573 


± 





069 






R c 


20.88 


± 





12 
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20.650 


± 
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R T 


20.84 


± 





13 






A 0,e 
^B 


0.0170 


± 





0049 






a0,m 

^B 


0.0164 


± 
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A 0,r 
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0.0240 


± 
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mf 


-92 


[GeV] 


91.1883 


± 





0084 


mf 
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[GeV] 
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± 





0043 


mf 


[GeV] 


91.1894 


± 





0038 
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correlation matrix 
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— U 
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U 
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[GeV] 
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Z 
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± 
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The common energy errors are shown in Table The agreement of the energy- related 
errors for the eleven parameter results is acceptable among the experiments. Noticeable are 
the differences in the covariance matrix elements between the Z mass values and a^', they 
were traced back to slightly different strategies in choosing a value of mz in the different 
time periods for calculations involving the hadronic pole cross section. 
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r z 


< 


Re 


ml " 92 


m 93 " 94 


95 

m z 


Aleph 


m z u ~ 


-92 


-0.0011 


-0.002 


0.0020 


0.0058 










-94 


-0.0007 


-0.003 


0.0016 


0.0006 


0.0028 






95 

m 2 




0.0002 


-0.002 


0.0003 


0.0003 


0.0013 


0.0015 


Delphi 




-92 


-0.0009 


0.0014 


0.0027 


0.0053 








mf- 


-94 


-0.0006 


-0.0020 


0.0014 


0.0008 


0.0027 






95 

m| 




-0.0002 


-0.0018 


0.0007 


0.0004 


0.0012 


0.0014 


L3 


m y z °- 


-92 


-0.0013 


0.0045 


0.0046 


0.0050 








m| 3 ~ 


-94 


-0.0007 


-0.0022 


-0.0007 


-0.0008 


0.0028 






95 

m z 




-0.0001 


-0.0018 


-0.0006 


-0.0009 


0.0012 


0.0015 


Opal 




-92 


-0.0011 


0.0025 


0.0036 


0.0058 








m 93 " 


-94 


-0.0010 


0.0034 


0.0034 


0.0006 


0.0030 






95 




0.0000 


-0.0014 


0.0013 


-0.0003 


0.0012 


0.0015 


average 




-92 


-0.0011 


0.0016 


0.0032 


0.0055 








m z 3 ~ 


-94 


-0.0008 


-0.0010 


0.0014 


0.0003 


0.0028 






95 




0.0000 


-0.0018 


-0.0004 


-0.0001 


0.0012 


0.0015 



Table 19: Common energy errors for multiple-mz fit, from individual data sets and average. The 
matrix elements for Tz, o~9 and R e and for the asymmetries are very similar to the nine-parameter 
case and can be taken from Table\l^ and\lo]. 
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